Nuclear Pairing in the T=0 channel revisited 
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Recent published data on the isoscalar gap in symmetric nuclear matter using the Paris force 
and the corresponding BHF single particle dispersion are corrected leading to an extremely high 
proton- neutron gap of A ~ 8 MeV at p ~ 0.5po- Arguments whether this value can be reduced due 
to screening effects are discussed. A density dependent delta interaction with cut off is adjusted so 
as to approximately reproduce the nuclear matter values with the Paris force. 



In a recent publication the possibility to reproduce 
the gap in nuclear matter, as obtained e.g. from the Paris 
NN force, by an effective density dependent zero range 
force, was investigated. Supplied with an energy cut off 
such effective forces turned indeed out to be able to re- 
produce very reasonably the gap values in the isospin 
T=0 and T=l channels over the whole relevant range of 
densities. The adjustments were performed on previously 
published solutions of the gap equation using Brueckner- 
Hartree-Fock results for the single particle spectra [^. 
Such effective forces may possess some analogies with 
similar ones frequently used in recent structure calcu- 
lation of superfluid nuclei Unfortunately, due to the 
subtleties connected with the numerical solution of the 
gap equation, the published results in the T=:0 channel 
were not accurate enough so that the corresponding gap 
is underestimated in by about 20%. It is the purpose 
of this note to give the corrected results for the gap in 
the T=0 channel and also to readjust the corresponding 
density-dependent (5-force. We also discuss again the is- 
sue whether screening affects the T=l and T=0 channels 
differently. 

In Fig. 1 we show the correct result for the isoscalar 
gap as obtained with the Paris force [Q using two in- 
dependent numerical codes. We also checked that the 
Argonne V14 force |^ gives practically the same result. 
What is striking is the giant gap value of ^ 8 MeV at 
maximum, which is of the same order as the Fermi energy 
at the corresponding density. Even around saturation, A 
is still of the order of several MeV. This is clearly a strong 
coupling situation as expected from the fact that at low 
density the n-p Cooper pair turns into the deuteron wave 
function Q. The above values are actually much more 
compatible with earlier calculations of the critical tem- 
perature in Ref. than the previous results 1 1] . Indeed, 
considering the usual relation A = 1.76 Tc , quantita- 
tive agreement between the results of and the ones in 
Fig. 1 is obtained. In order to obtain an estimate of the 



typical magnitude of the isoscalar gap in a finite nucleus, 
we apply the local density approximation and average the 
local gap over the density at the Fermi energy. This pro- 
cedure has given reliable estimates of the average energy 
dependent gap in the isovector channel . We therefore 
calculate 
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FIG. 1. Pairing gap versus Fermi momentum for symmetric 
nuclear matter in the T=0 channel from the Paris potential. 
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where the local Fermi momentum is defined as 
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with /i the chemical potential. We take the same single 
particle potential V{r) as in Q and the result for e.g. 
N=:Z=35 is that A is of the order of 3 MeV. Compared 
to the neutron-neutron and proton-proton channels this 
is a very high value. 

We already discussed in and show again in Fig. 2 
that the use of the Paris force in conjunction with the 
k-mass, m* /m, yields gap values as a function of density 
which are globally very similar to the ones of the Gogny 
force for T=l and therefore, the use of a bare force seems 
not unreasonable in the T=l channel. The fact that A 
for T=l drops off quite a bit faster close to saturation 
for the Paris force than for the Gogny DIS force may be 
attenuated in a finite nucleus to quite some extent, since 
a certain averaging over all densities p < po takes place. 
Therefore the needed medium renormalization of the bare 
force seems to be of minor importance in the T=l chan- 
nel^. However, the situation may not be the same for 
T=0 pairing. The extremely strong T=0 pairing stems 
essentially from the fact that with respect to the T=l 
channel the tensor force is acting additionally. Without 
the tensor force np (T=:0) and nn (T=l) pairing would 
be of comparable magnitude. The screening of the tensor 
force in the medium is, however, still a controversial sub- 
ject 1^. On the other hand, even for very low densities 
where screening should not be so important, T=0 pairing 
remains strong. Therefore, there may be a good chance 
that the new heavier exotic nuclei with N — Z experi- 
ence quite pronounced np superfluidity. This may well 
be the cause for the so called Wigner energy of the nu- 
clear mass formula, since it can be shown pO| that away 
from symmetric nuclei, T=0 pairing looses very quickly 
its strength. 

Let us now proceed to the readjustment of the effective 
T=0 delta force. We use the standard ansatz lll],ni| 
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With the above density-dependent zero range force, the 
gap equation reads 
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FIG. 2. Pairing gap Af in tlie ^So channel in symmetric 
nuclear matter calculated with the Gogny force DIS compared 
with results from the Paris force. 
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In Fig. 3 we present two fits for the above ansatz, 
one of the fits is obtained from the following parameters: 
a = 0.2, 77 = —0.10 and a cut off energy ec = 60 MeV 
(see Ref. using the effective mass m* /m as obtained 
from the Gogny force. The other fit is obtained by us- 
ing a bare mass and parameters a — 0.90, rj = 0.40 and 
ec = 60 MeV. 

As one can see in Fig. 3, the fit obtained using the 
bare mass is able to reproduce the microscopic calcula- 
tion up to the highest values of kp [kp ~ 1.7 fm~^), while 
the fit obtained using the effective mass breaks down at 
lower densities corresponding to kp ^ 1.35 fm~^. The 
reason for this different behavior can be traced back to 
the dependence on the effective mass inside the integral 
of the gap equation. It turns out that in order to get 
a solution of the gap equation (^, the energy cut off 
ec should be larger than the Fermi energy ep. Other- 
wise no value of A satisfies the equation. In the case 
of the energy cut off used in Figs. 3 and 4 (ec = 60 



'^Of course it cannot be excluded that the medium com- 
pletely re-shufHes the distribution of gap values, still repro- 
ducing experimental pairing phenomena in finite nuclei 
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MeV), the largest kp reachable is kp ~ 1.7 fm~^ when 
bare masses are used, but only kp ~ 1.35 when effective 
masses are used instead. Therefore, we plot in Figs. 3 
and 4, the fits obtained only up to those values of kp, 
when m*/m ^ 1. Nevertheless, the fits cover all the 
physically relevant range of densities form zero to satu- 
ration (po — 0.16 fm~^). 
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FIG. 3. T=0 pairing gap in nuclear matter. The dots are 
the results obtained for the Paris potential. The curves are fits 
with Eq.(^ using an energy cut off ec ~ 60 MeV, i;o = —480 
MeV fm^, and different parameters for the fit with effective 
mass m* (solid line, rj — —0.10, a — 0.20) and for the fit with 
the bare mass (dashed line, rj = 0.40, a — 0.90). 



In principle, in the T=0 channel, vq should be cho- 
sen such that the deuteron binding energy is reproduced 
in free space. However, we have found that with this 
condition the fit obtained is very poor. Therefore, for 
a given energy cut ec, we vary the parameter vq from 
the value that produces a bound state at zero energy 
vo = —(h'^/m)(2n'^/\/2rnec), up to the value that pro- 
duces the bound state at the deuteron energy and 
choose the best fit. The fits in Fig. 3 have been obtained 
with Vq = —480 MeV fm"^ as it corresponds to a bound 
state at zero energy. This reduces the value of the gap at 
low densities but improves significantly the fit at higher 
energies. On the other hand, as we shall see in Fig. 4, 
the value of ?;o is chosen between the two extreme val- 
ues considered, bound state at zero energy and at the 
deuteron energy. In any case, the values used for vq are 
quoted in each case. 
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FIG. 4. Same as in Fig. 3 but suppressing the density de- 
pendence (?) = 0) and using vo = —530 MeV fm^. 

In Fig. 4 we present a similar fit for the case with 
m* /m 1, however, suppressing the density dependence 
completely, that is = 0. Since this parameter was al- 
ready small for the case in Fig. 3 the fit is still acceptable 
and only a slight deterioration at the low density end is 
visible. Let us mention also that the use of the bare mass 
m* — m allows an excellent fit of the microscopically cal- 
culated gap values at all densities (see Fig. 3). However, 
realistic calculations of finite nuclei are rarely performed 
with the bare nucleon mass. 

As a first guess we may try to use the effective pairing 
force obtained with the present fit also for finite nucleus 
calculations. This will give a rough account of whether 
the use of a bare force in a finite nucleus is at all rea- 
sonable in the T=0 channel. We would, however, like 
to point out that the expression of Eq.(||) for finite nu- 
clei may not give precise reproduction of the results one 
would obtain with a direct use of the Paris force in the 
gap equation. Indeed, in the mean time, we compared in 
the T=l channel the results of the genuine Gogny force 
and its density dependent (5-force substitute elaborated 
in 111 in a half infinite matter calculation [Q. Prelim- 
inary results show that the detailed surface dependence 
of the gap and of the anomalous density seem to be quite 
different in both cases. However, integrated quantities 
like the correlation energy may still be rather similar. 

Of course, it should be interesting for the future to de- 
rive also an effective finite range force in the T=0 channel 
which is as efficient as the Gogny force for T=l pairing. 
In fact the Gogny force has never been used for np pair- 
ing. However, since in this channel the density dependent 
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zero range force enters, one has to introduce an additional 
cut off which is an unknown adjustable parameter. 

In summary we give corrected values of the np (T=0) 
gap in nuclear matter using the Paris force together with 
Brueckner-Hartree-Fock single-particle energies. An ex- 
tremely high value of A 8 MeV at p ^ 0.5po is ob- 
tained, leading to a gap value in finite nuclei of ~ 3 
MeV. Argmiients are advanced that the pairing force in 
the T=0 channel may be more strongly screened than in 
the T=l channel. We then adjust a density-dependent 
5— force to the nuclear matter gap values. The fit is rea- 
sonably successful for densities below saturation. 
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